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Abstract 

We present a model where an unstable MeV Majorana tan neutrino can naturally 
reconcile the cold dark matter model (CDM) with cosmological observations of large 
and small scale density fluctuations and, simultaneously, with data on solar and atmo- 
spheric neutrinos. The solar neutrino deficit is explained through long wavelength, so- 
called just-so oscillations involving conversions of Ve into both and a sterile species 
I's , while atmospheric neutrino data are explained through to fg conversions. Fu- 
ture long baseline neutrino oscillation experiments, as well as some reactor experiments 
will test this hypothesis. The model is based on the spontaneous violation of a global 
lepton number symmetry at the weak scale. This symmetry plays a key role in generat- 
ing the cosmologically required decay of the Ur with lifetime t^^ ~ 10^ — 10^ seconds, as 
well as the masses and oscillations of the three light neutrinos I'e , t'/i and I's required in 
order to account for solar and atmospheric neutrino data. It also leads to the invisibly 
decaying higgs signature that can be searched at LEP and future particle colliders. 
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1 Introduction 



A tau neutrino with a mass in the MeV range is an interesting possibihty to consider for 
two different reasons. On experimental side such a neutrino is within the range of the 
detectabihty, for example at a tau-charm factory [Q, |^. On the other hand, if such neutrino 
decays before the matter dominance epoch, its decay products could then add energy to 
the radiation thereby delaying the time at which the matter and radiation contributions to 
the energy density of the universe become equal. Such delay would allow one to reduce the 
density fluctuations at the smaller scales purely within the standard cold dark matter 
scenario @], and could reconcile the large scale fluctuations observed by COBE with the 
earlier observations such as those of IRAS on the fluctuations at smaller scales. An MeV 
v-r may, however, conflict with the big-bang nucleosynthesis picture [0]. This conflict can be 
avoided in two ways. 

If the tau neutrino has a strong coupling to a light spin zero Goldstone boson - called 
majoron and denoted J - with a typical strength gj > 10~^ , then the annihilation of Ur pairs 
to majorons could reduce their number density sufficiently so as to be consistent with the 
nucleosynthesis bound [Q. In spite of this reduction, subsequent decays could increase the 
energy density of the radiation enough to reconcile COBE with IRAS data. This possibility 
has been recently investigated in the context of a specific doublet majoron model 
where an upper bound on the required life of 10^ seconds was set. 

Alternatively, even if the majoron coupling to the tau neutrino is not so strong it may 
be possible to reconcile the nucleosynthesis constraints with the MeV Ur hypothesis if its 
decay involves the electron neutrino. The decay z/g could be captured by neutrons so as to 
reduce the resulting yield of primordial helium. In this case, a u^- with mass of a few MeV 



and lifetime in the range 10-50 seconds has been advocated in ref. ^0 



An MeV tau neutrino, though cosmologically interesting, does not obviously fit with the 
data on solar and atmospheric neutrinos [|rT[], for which neutrino oscillations with quite small 
mass differences provide the most plausible solutions |1T2| , p!3| . However, so far all attempts 



14 1 to obtain an MeV Ur with the lifetime required to revive the cold dark matter picture 



have ignored solar as well as atmospheric neutrino data. It is desirable to develop a coherent 
model which not only fits COBE and IRAS observations, but also provides solutions to the 
solar and atmospheric neutrino problems. 

In this note we realize the Mev tau neutrino hypothesis in a model that can naturally 
reconcile the cosmological data on primordial density fluctuations with an explanation of the 
solar and atmospheric neutrino deficits through neutrino oscillations. The simplest way to 
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do this is to assume the few MeV Vr to be a majorana neutrino and not a Dirac neutrino as 
assumed in [|10]. Indeed, an Mev Vr can not pair up with or with z/g in order to form a 
Dirac state, because of the laboratory bounds on the masses of z/^ or z/g. As a result, an 
MeV Dirac state would be obtained only by pairing off the two-component z/^ with a sterile 
neutrino state of the same mass. In this case at least three other light neutrino species, 
one of which should also be sterile, would be required in order to fit together the data on 
solar and atmospheric neutrino oscillations. This follows from the fact that in this case the 
oscillations of z/^ or z/g into the MeV z/^ cannot solve the atmospheric or the solar neutrino 
problems and, on the other hand, the Ve — ^'/i oscillations could solve either but not both, 
since they require quite different values for the corresponding (mass)^ differences JT^ . 



Hence, the most economical way to reconcile COBE and IRAS observations with the 
cold dark matter picture and with solar and atmospheric neutrino data should involve the 
presence of just one very light sterile neutrino and just one two-component MeV state: the 
majorana tau neutrino. 

To construct a consistent and appealing theoretical model is a non trivial task. Apart 
from naturally relating the required smallness of the sterile neutrino mass to a suitable 
symmetry, one has to obey a number of constraints: 

(i) The mass and mixing of the very light sterile neutrino i>s should not conflict with the 



constraints coming from the nucleosynthesis [|T5l and supernovae |]16 



{ii) The model should contain three different mass scales namely m,^^ ~ MeV to account 
for the tau neutrino mass, Aq ~ 10^^ eV ^ or 10~^° eV ^ in order to account for the solar 



neutrino deficit and ~ 10^^ eV ^ to solve the atmospheric neutrino problem. 

(iii) The u^- should decay with lifetime in the range 10^ — 10^ seconds [§] or 10 — 50 seconds 
p!0| . In the former case, it should couple strongly to majoron, while in the latter case its 



decay products should produce z/g. 

(iv) The couplings of the majoron should be strong enough to satisfy [iii), but this coupling 
should not result in excessive energy loss through majoron emission inside a supernova [|T^. 

We present below a model which successfully meets all these requirements and discuss 
its most salient features. 
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2 The model 



Our model is based on the triplet plus singlet majoron scheme |jT8| and contains three right 
handed neutrinos, one of which is kept light by the imposed global symmetry. This way of 
keeping the sterile neutrino light has been already used in variety of models which tried to 
accommodate the possible existence of a 17 keV neutrino state [jl9| or which try to solve the 
solar, atmospheric and the dark matter problems simultaneously ||12| , pO| |. 



We replace the lepton number symmetry of the original singlet majoron model by a 
generation dependent global symmetry U{1)g under which the various fields transform in 
the manner shown in table 1. The generation dependent symmetry serves two purposes. It 
keeps the sterile light and it leads to the decay of the tau neutrino into lighter neutrinos 
plus a majoron The SU{2) x U{1) x U{1)g invariant couplings of the neutrinos are 

given by: 



< 

< a > 2 



+ H.C. 



The above Yukawa interactions lead to the following neutrino masses: 

1 



-ul' CM.ul + H.C. 



(2) 



where Mj^ is a 6 x 6 matrix having the following form in the left-handed basis v = 



fz/^, i/e, u^, z/^, z/;^, u^)^: 
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We assume the bare mass M to be much greater than other scales appearing in eq. @, 
as in the seesaw model. These two heavy states can be diagonalized out leading to four light 
^In principle, these problems can be solved without invoking a light sterile state if neutrinos are almost 



degenerate ||20| , 21 1. However, in our present case one is obliged to introduce a light sterile neutrino if the 
mass of i/r is to lie in the MeV range. 
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states, one of which is massless, from the form of eq. @. The effective 4x4 neutrino mass 
matrix obtained in the seesaw approximation takes the following form: 



m 



(3cose I3sme 



13 cose 
/3sin^ 



a cos 6 
a sin 9 



(4) 



where, 



and 



^ a cos 6* a sin 6' rrij,^ y 



n2 — / \2t 2 I 2\ 2 _ /^3n,2/ 2 , 2\ 



■M' 



tan 6' 



m2 
mi 



(5) 



(6) 



Note that the v^r is not allowed to receive a large mass and is kept light after the seesaw 
mechanism. As already mentioned, one of the four light states described by w^e// is in fact 
massless, Ai = 0. The other three are massive with eigenvalues A2,3,4 approximately given 
by: 

„2 



A, 



a 



2m,. 



(7) 



A3 ^ /5 + 
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X4 



These eigenvalues nicely reproduce the hierarchical scales required for our purposes. Because 
of the chosen quantum numbers with respect to U{1)g, the Uj. is the only state to receive the 
mass from the SU (2) triplet fields and could be in the MeV range. The parameters a and 
/3 could be much smaller than m,^^ if the seesaw mass scale M is chosen appropriately large. 
In this case, the matrix me// itself has a seesaw structure. The effective matrix describing 
the mixing of z/g,^ and the sterile neutrino (z/^ = UeR) is easily seen to posses an approximate 
Lg — — L(. symmetry. This leads to a pair of almost degenerate states A2,3. Their mass 
provides the (mass)^ difference 



A, 



(9) 



The approximate symmetry and hence degeneracy among two of the neutrinos is broken 
by the corrections O {:^) arising out of the vacuum expectation value (VEV) of the a 
field which breaks the global symmetry and generates the majoron. These corrections are 
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naturally small because of the second seesaw mechanism and they provide another (mass)^ 
difference 



A,. = A - A 



m,. 



(10) 



For the values a ~ /? <^ rriiy^ one naturally obtains the hierarchy ^ A^. 
The mixing among the four light states is specified by the matrix: 



U rUeff = diag.{0, 171,^2, m„^,m„^) 



fill 



where m^- = |Aj| and U is defined by 

/ 





1 

V2 



V2 



13a 



— sin 61 ^ cos 9 — cos 9 



■ cos( 



v 



COS 9 sin I 




V2 



sin 9 — — sin 9 



The weak eigenstates are related to the mass eigenstates as follows: 

1 / .X /3a 
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(14) 



This gives rise to the following pattern for neutrino oscillations. Consider first the limit in 
which As is neglected in comparison to A^. In this limit, there are no oscillations involving 
the sterile state us . On the other hand the z/g and z/^ oscillate among themselves with the 
probability: 

P^, = sm' 29 sin^ (^) (15) 

These oscillations could give rise to an explanation of the observed depletion in the atmo- 
spheric neutrino flux if the parameters lie in the range [0, |T^ 



sin^ 29 = 0.35 - 0.8 ; A^ = (0.3 - 2) x 10"^ eV^ 



(16) 



When the smaller (mass)^ difference As is turned on, the and the start oscillating 
into sterile state i>s . The probabilities for the z/g oscillations averaged over the shorter 
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atmospheric neutrino oscillation scale set by is given by 



Pe,{t) = ^(3 + COs(||; 

Pes{t) = ^fl-COsl^)) (17) 



2 V "2E 



Since the mixing angle involving z/2 and z/3 is 45°, the nonadiabatic MSW solution |2^ to the 
solar neutrino problem is not feasible in the present case. However, the relevant mass scale 
As could be naturally very small so that the solar neutrino flux may get depleted through 



the long wavelength vacuum oscillations p4|. For example, m^^ = 5 MeV, /? = 0.1 eV and 
a = 0.05 eV would give As = 10~^°eV^^ and = 10~^eV^. One sees that similar values 
of a and (3 naturally result in a very large hierarchy between the solar and atmospheric 
scales. On the other hand, if one chooses the Dirac masses mi 2,3 in the GeV range, then the 
required values of a, (3 can be obtained by choosing the bare mass M for the right handed 
neutrino around the intermediate scale 10^^ GeV. One concludes that the required values of 
A^ and A^ do follow for natural choices of the parameters. 

Note that, although it has been shown that Ve to z/5 oscillations where z/5 is a sterile 
state are disfavored by the combined Homestake and Kamiokande data |]2S|, in the present 
case the situation is more complex, since not only the sterile state z/5 but also are involved. 
Since both relevant mixing angles are large |^ a phenomenologically consistent solution should 
exist. In order to determine its parameters more sharply a more detailed analysis of the 
existing solar neutrino data for the present case where both and vs take part in the solar 
neutrino oscillations would certainly be desirable. 



3 Cosmology 

Let us now turn to the cosmological aspects of the model. A stringent constraint on this 
scenario comes from primordial big bang nucleosynthesis, which requires that the effective 
number of degrees of freedom g^ffiT) contributing to the energy density of the universe at 
the nucleosynthesis time (~ 1 second) be less than 11.3 0. Here QeffiT) has been defined 
in terms of the total energy density as p = ^geffiT)T'^ which includes the contribution 
§A possible restriction on large angle neutrino mixing has been argued to follow from the observed energy 



spectra of i^e from supernova SN1987A |26 . However, ref. considered only the simplest case of two flavour 
to mixing, whereas in the present case one has three neutrino species one of which is sterile, so that 
those arguments do not directly apply. 
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of the relativistic species as well as that of the nonrelativistic Ur- The latter could violate 
this bound if it had only the conventional weak interactions. However, the presence of the 
majoron in our model alters the situation. The majoron in the model is easily seen to be 

j^,,.^T,+ '-^H (18) 

u uv 

where u,v, and to denote {\/2 times) the vacuum expectation values of the a, 0° and T° 
scalar multiplets, respectively, while the suffix / denotes the corresponding imaginary parts. 
Note that because of the hierarchy uj ~ MeV ^ m ~ lOOGeV^, the invisible decay of Z 
to the majoron is enormously suppressed, in accordance with LEP data [^, unlike in the 



purely triplet majoron scheme. The Vr couples to the majoron dominantly through its triplet 
admixture. Using eq. (|lD and eq. (^) this coupling may be given as follows 

{Tfi (y. V 1 I 

^^Z/fiCz/4L + Z/J^C(Z/2L - Wzl) + {l^2L " I^ZL^CVaL \ + H.C. (19) 

u V2u L J J 

The Vr ~ has a diagonal coupling to the majoron given by gj = For sufficiently 

large gj the contribution of at the time of nucleosynthesis can be suppressed through 
i/t annihilation. In order for this to happen, one requires 

Teqi = ^ m,^Y < 0.13 MeV (20) 

which corresponds to geff{T ~ IMeV) < 11.3. Here Y determines the abundance of the 
nonrelativistic Ut-. Using the standard expression for Y we obtain 



« (2.66 X [-^) {^^) (f ) (21) 

where 

X/ = In^ - (n + l/2)lnlnz (22) 

and 

z = 0.038(n + l)g/gl/^Mpim,^(To (23) 
where n = 1, Mpi = 1.2 x 10^^ GeV is the Planck mass, g = 2 and, in our case [@], (?* = 10. 

2 

Finally, do = |^ denotes the Ut- annihilation cross section. 

Physically, Teqi represents the temperature at which the contribution of the Ur starts 
dominating over that of radiation. As would be expected, this should happen after the 
nucleosynthesis epoch. After Teqi, the u^. would dominate the energy density of the universe 
until it decays. This decay would make the universe radiation dominated and this could serve 
to delay the time at which matter starts dominating again. If one fixes the overall scale of the 
power spectrum of the density fluctuations in the cold dark matter scenario from the large 
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scale COBE observations, then one can also fit the small scale density fluctuations provided 
one chooses geff{TEQ2) ~ ^ ^ 3.36 ||, jT^l- Here geff(TEQ2) determines the total contribution 
in radiation at the temperature Teq2 when the radiation and matter densities again become 
equal. The i^eji^ti, majoron and photons contribute the amount 3.17 to geff{TEQ2)- The 
remaining must come from the decay products of the Ur- This can be translated |0] into a 
constraint on the lifetime: 

/ 7^0 \2 

r,, ^ 7.2 X 10^9 —- sec (24) 

where T° = 2.73 is the present temperature of the universe. Since Teqi < 0.13 MeV from 
nucleosynthesis, r^^ > 2 x 10^ sec. On the other hand, Teqi should be higher than the 
temp Teq2 ~ few eV of the conventional matter radiation equahty, giving an upper bound 
< 10^^ sec. In fact, requiring that the Ur contribution to the present total energy density 
of the universe is not too large gives a more stringent upper bound on r^^ as a function of 
the mass m^^ |2^. For example, for m^^ ~ 1 MeV, T,y^ < 10^ sec. 

The above considerations show that a T,y^ in the range ~ 10^ — 10® seconds would be 
able to delay matter radiation equality without conflicting with the nucleosynthesis picture. 

We now show how the lifetime required by cosmology can naturally occur in the present 
model. For this note that eq. (|19D gives rise to the decay of Ur to the lighter states 1^2,3 plus 
a majoron, with decay width given by 

2 



a 



47r \u, 

Using the expressions for and this can be written as 



(25) 



= (1.6 X lOW) (^) [^^) (^) ] (26) 

It follows from the above equation that one can simultaneously accommodate the solar and 
atmospheric neutrino deflcits and also have a tau neutrino decay with the required lifetime 
if the singlet VEV u is chosen around 50 — 200 GeV. In Fig.l we display the region of 
mass and the singlet VEV u allowed by the various constraints. The solid curves a 
and b correspond to (^e// = 10.86 and 11.3, respectively. The regions to the right of these 
curves would be disallowed by nucleosynthesis. For illustration, we give curves c and d 
corresponding to v.^ lifetimes 10^ and 10^ seconds, respectively. 

Another constraint on the model comes from the supernova. The MeV tau neutrino 
with relatively large couplings to the majoron may cause supernova to rapidly loose energy 



through majoron emission. The relevant constraints have been looked at in detail in ref. [|T7 
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in the case of the simplest singlet majoron model. In our case the relevant coupling is the 
first in eq. ([T9|) . Since, this coupling is similar (apart from a factor 2) to that of the singlet 
majoron model we can adopt the analysis of ref. [jl^. In our case, the following values of u 
and m^^ are seen to be disallowed 

5.7xl0"^< < 0.82 27 

It follows from the above two equations that it is indeed possible to obtain the Ur lifetime 
in the required range and solve the solar and atmospheric neutrino problem for parameter 
choices lying outside the range forbidden by the supernova. 

The time structure of the SN87A antineutrino pulse may also be used to constrain the 
lifetime of z/^. The decay of a massive emitted in the supernova on the way could lead 



to a delayed signal in the detectors. The absence of such a signal has been used in [30| and 



subsequently in [31] to put an upper bound on the lifetime. A rough estimate of the 
resulting lifetime has been given in ref. to be ~ 300 seconds for an Mev z/^ decaying 
by majoron emission. While this is consistent with the bound we have obtained of 2 x 10^ 
seconds, there is considerable room to relax this supernova bound. Firstly, only a fraction 
cos^ 9 of Ur decay to Ue in the present case. Moreover, original z/^ flux may also be suppressed 
by the Boltzman factor for larger masses m,y^ > lOMeV. These could result in considerable 
weakening of the upper bound. Finally, in the derivation of the supernova bound one would 
have to be careful in differentiating between 1/^ and z/^ decays to z/g. 

The presence of a light sterile neutrino mixing with ordinary neutrinos is in general 
constrained by nucleosynthesis |T^. Neutrino oscillations in the early universe would bring 
the sterile neutrinos into equilibrium at the time of nucleosynthesis with the active ones. 
The relevant oscillation scale in the present case is ~ 10~^°eV^. The corresponding 
wavelength is too large to equilibrate the sterile neutrinos in the early universe. In fact the 
constraints on the relevant mass scale is trivially satisfied in our case. 



4 Discussion 



We have attempted in this paper to provide a coherent explanation of quite distinct phe- 
nomena in neutrino physics. The example given here is able to resurrect the cold dark 
matter picture of structure formation by making it consistent with both COBE and IRAS 
observations of primordial desnsity fluctuations through a Ur of few MeV mass and lifetime 
in the range 10^ — 10^ sec. It also explains the solar and atmospheric neutrino problem in 
terms of mixing among three light neutrinos z/g , and us . The scheme presented here 
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is not unique but is certainly most economical from the point of view of explaining various 
phenomena mentioned above. More importantly, it contains predictions which can be tested 
in the near future in long baseline neutrino oscillation experiments, as well as some reactor 
neutrino experiments. Moreover, due to the presence of the weak scale majoron, the present 
model allows for a distinctive signature of the invisibly decaying Higgs boson h J J 
which could substantially affect higgs boson search strategies at LEP NLC [|4l) as well 



as LHC |5|. 

Finally, notice that, since is much heavier than other neutrinos and has very small 
mixing with them [O [a/niy^ ~ 10^'')) there are no oscillations involving the and exper- 
iments such as CHORUS and NOMAD looking for the z/^ — Vj. oscillations should show 
negative result. This feature is a generic expectation in any model having in the MeV 
range. It also implies a strong suppression in the neutrinoless double beta decay rate. 
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for fruitful correspondence. 



10 





SU{2) 


Y 


G 




2 


-1 


-1 




2 


-1 


1 




1 


-2 


-1 


tr 


1 


-2 


1 


l^eR 


1 





-2 


I^HR 


1 





-1 


^tR 


1 





1 





2 


-1 





T 


3 


2 


-2 


(7 








1 



Table 1: SU{2) ® U{1)y ® U{1)g assignments of the leptons and Higgs scalars. Quarks are 
U{1)g singlets. 

Figure Caption 

Fig.l shows the region of u^- mass and hfetimes as a function of the singlet VEV u which 
are allowed by the various constraints. The solid curves a and b illustrate the nucleosynthesis 
contraints corresponding to geff = 10.86 and 11.3, respectively. The regions to the right of 
these curves would be forbidden. Curves c and d correspond to i/^ lifetimes 10^ and 10^ 
seconds, respectively. 
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